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The detection of neutrinos from SN 1987A by the Kamiokande-II and Irvine-Michigan-Brookhaven detec- 
tors provided the first glimpse of core collapse in a supernova, complementing the optical observations and 
confirming our basic understanding of the mechanism behind the explosion. One long-standing puzzle is that, 
when fitted with thermal spectra, the two independent detections do not seem to agree with either each other or 
typical theoretical expectations. We assess the compatibility of the two data sets in a model-independent way 
and show that they can be reconciled if one avoids any bias on the neutrino spectrum stemming from theoretical 
conjecture. We reconstruct the neutrino spectrum from SN 1987 A directly from the data through nonparametric 
inferential statistical methods and present predictions for the Diffuse Supernova Neutrino Background based on 
SN 1987A data. We show that this prediction cannot be too small (especially in the 10-18 MeV range), since 
the majority of the detected events from SN 1987A were above 18 MeV (including 6 events above 35 MeV), 
suggesting an imminent detection in operational and planned detectors. 

PACS numbers: 97.60.Bw, 98.70. Vc, 95.85.Ry, 14.60.Pq 



I. INTRODUCTION 

Core-collapse supernova explosions are how massive stars 
end their lives JH. Once their nuclear fuel is exhausted, they 
produce a neutron star or black hole, recycle their ashes into 
cosmic dust to seed subsequent starbursts and spread then- 
heavy elements 10, [HQ], and flood the universe with a burst 
of neutrinos. Supernova explosions are very prolific neutrino 
producers, since the enormous amount of gravitational bind- 
ing energy, liberated when the core of a massive star col- 
lapses, can effectively be disseminated only by neutrinos. In 
simulations, the rebounding collapse of the core often pro- 
duces a shock that stalls, which may indicate the necessity 
of shock revival by ener gy t ransfer from the intense neutrino 
flux Sll, SHIS [Unl El. While neutrinos can be cru- 
cial for diagnosing a successful explosion and revealing the 
conditions of the core, their elusive nature makes their detec- 
tion challenging. 

The detection of neutrinos from SN 1987A in the Large 
Magellanic Cloud (at a distance of D ~ 50 kpc) by 
the Kamiokande-II (Kam-II) IU4l [1511 and Irvine-Michigan- 
Brookhaven (1MB) flrjl [nil detectors transformed super- 
nova physics from purely hypothesis- and simulation-driven 
to discovery-driven science. One puzzling feature of the 
SN 1987 A data is that the neutrinos detected by the 1MB 
detector were seemingly more energetic than those detected 
by the Kam-II detector, which were clustered at low energies. 
Despite the uncertainties associated with the low statistics and 
the differences of the detector properties (like efficiency and 
size), it seems hard to accommodate the two detections with a 
common quasi-thermal neutrino spectrum lfl8ll . This incom- 
patibility between the data sets and the theory can be allevi- 
ated by assuming that at least one experiment had an anoma- 
lous statistical fluctuation or underestimated systematics lfl9Tl . 
Here, we assume that both experiments observed statistically 
probable outcomes, as they were exposed to the same spec- 
trum, and accept the reported events, uncertainties and overall 
efficiencies at face value. 



The principal focus of our study is to assess the compatibil- 
ity of the Kam-II and 1MB data sets in a model-independent 
fashion, and to study the immediate implications for neutrino 
detection from Galactic and cosmic supernovae. We begin 
with a brief formulation of neutrino detection from a nearby 
supernova and summarize the main ingredients. Next, we 
determine the relative prospects of detecting neutrinos from 
SN 1987A in either Kam-II or 1MB, based on their reported 
energy-dependent efficiencies and the numbers of targets in 
their fiducial volumes. We compare this relative probability 
of detection to the actual detected positrons for Kam-II and 
1MB and conclude that there is no significant discrepancy be- 
tween the data sets unless an a priori spectral shape is forced 
to fit the two data sets simultaneously. We establish that the 
Kam-II and 1MB detectors mainly probed different energy do- 
mains, the former testing the low energy part of the spectrum 
and the latter testing the high energy part, and, in fact, only 
a proper combination of Kam-II or 1MB data sets provides a 
sensitive probe over all energy ranges. 

Using nonparametric (distribution-free) inferential statisti- 
cal methods [20], we determine the structure of the underlying 
spectrum directly from the data. Since these methods make 
no a priori assumptions and do not rely on parameter estima- 
tion, they allow for more efficient processing of small data 
samples. Then, we present the incoming neutrino spectrum 
inferred directly from the data, which would be an effective 
spectrum received on Earth after effects modifying the over- 
all shape are taken into account. Finally, we provide predic- 
tions for the Diffuse Supernova Neutrino Background based 
on the SN 1987A data (as Refs. fH HI have also recently 
done through more conventional methods), which turn out to 
be nearly as large as typical predictions in the literature based 
on supernova models. Our study follows a transparent ap- 
proach to extract the targeted information directly from the 
SN 1987 A data, skipping the intermediate stages of multipa- 
rameter function fitting, commonly used in the past literature. 
The methods we present avoid over-interpretation of limited 
data, while still allowing for a substantial analysis. 
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II. DETECTION OF SUPERNOVA NEUTRINOS 

We concentrate on the v e flux received on Earth from a su- 
pernova since the detectors are particularly sensitive to the in- 
verse beta decay, P e + p — » n + e + (we ignore possible small 
contributions due v + e~ or ^+ le O). We parametrize the v e 
number spectrum as (j>{E v ) = [L u /Eq] ip(E„), where <p{E„) 
is a normalized function (/ <p(E„) dE u = 1) describing the 
overall spectral shape with no prior assumptions. Here Eq 
is the average energy ((E„) = J E v <p(E„) dE v = Eq) and 
L v is the time-integrated luminosity (i.e., the total number of 
neutrinos is L v jEo) in this flavor. In the energy regime of 
interest, the cross section of inverse beta decay at the lowest 
orders in inverse nucleon mass, 1/M, is 

a{E v ) ~ 9.5xl(T 44 cm 2 (1 -6E„/M) [{E v - A)/MeV] 2 , 

(1) 

expressed in terms of the neutrino energy, E v B23I1 . The de- 
tected positron energy (visible energy), E + , can be related 
to the incoming neutrino energy as: E + = E v — A, where 
A = M n — M p is the nucleon mass difference (the higher 
order corrections are much smaller than uncertainties associ- 
ated with energy resolution). The positron spectrum in a given 
detector, as a function of visible energy for a supernova at a 
distance D, can be cast as 

*(E+) = S^-a(E u )ct>(E u ), (2) 

where the effective number of targets available in the detector 
is N(E + ) — N t r](E + ) as a function of visible (positron) en- 
ergy, in which N t stands for the number of free-proton targets 
weighted with the up-time fraction of the detector and r/(E+) 
is the detector efficiency function. 



III. COMPATIBILITY OF THE DATA SETS 

In the top panel of Fig. Q] we present the effective num- 
ber of targets available in the Kam-II and 1MB detectors by 
normalizing the targets in the fiducial volume of each detec- 
tor by the up-time and reported energy-dependent efficien- 
cies dHH, N(E+) = N t r)(E + ). The combined effective 
number of targets for both detectors together is also shown 
(in which they are treated as one big detector with energy- 
dependent efficiency, covering the full energy range sampled 
by the data). The number of free-proton targets in the fiducial 
volume of each detector was Nt 1.43 x 10 32 for Kam-II 
and N t =i 4.55 x 10 32 for 1MB [24]. Despite being — 3 times 
larger in size, 1MB typically had lower efficiency, especially 
at lower energies. Moreover, the 1MB detector was briefly of- 
fline whenever a muon passed through the detector, reducing 
the average up- time by ~13% lfl8tl . We emphasize that while 
all previous analyses did take the separate energy-dependent 
efficiencies of Kam-II and 1MB into account, our approach of 
considering an effective combined detector with a more com- 
plicated energy-dependent efficiency is new. 

Note that the efficiency of Kam-II (and most modern neu- 
trino detectors) is basically a step function over the detector 



threshold. Kam-II had a low threshold of E t hr ~ 7.5 MeV 
and its efficiency quickly rose, starting at ~ 5 MeV, to its 
asymptotic value. This is because the efficiency function 
around the threshold energy is non-trivial due to the proba- 
bilistic nature of the detection process (i.e., in principle, a low 
energy neutrino could trigger a greater number of photomul- 
tipliers and be detected). At the time of SN 1987A, almost 
a quarter of the photomultipliers in the 1MB detector were 
offline due to a power failure, resulting in a distinct energy 
dependence due to the geometrical modifications. Thus, 1MB 
had a higher threshold of E t hr ~ 19 MeV and rose gradually, 
starting at ~ 16 MeV. At low energies, the yields will be dom- 
inated by Kam-II, due to its much greater efficiency, while at 
high energies, the yields will be dominated by 1MB, due to its 
greater size. 
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FIG. 1: Top Panel: The effective number of targets (free protons) 
available in the Kam-II and 1MB detectors and their combination 
as a function of the visible (positron) energy E+ (i.e., protons in 
the fiducial volume of each detector normalized by the up-time and 
energy-dependent efficiencies). Middle Panel: The fraction of to- 
tal detected events expected to be assigned to either Kam-II or 1MB 
(independent of incoming neutrino spectrum) as a function of the 
visible energy. Bottom Panel: Positrons detected during SN 1987A 
by Kam-II (11 filled-diamonds) and 1MB (8 open-circles) with their 
recorded energies (Ei) and experimental uncertainties (AEi). 
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When both the efficiencies and numbers of target protons 
in the fiducial volumes of the two detectors are taken into ac- 
count simultaneously, it becomes apparent why Kam-II regis- 
ters all the positrons below 16 MeV, while 1MB is expected 
to dominate over 30 MeV. An informative perspective in com- 
paring the effectiveness of Kam-II and 1MB is presented in 
the middle panel of Fig. [T] This shows the fraction of total 
positrons expected to be recorded by either of the experiments 
(i.e., the ratio of effective number of targets in each detector 
to the combined number of targets) as a function of visible 
energy. Note that this is independent of the shape of the in- 
coming neutrino spectrum or the cross section, since it only 
compares the effectiveness of the detectors relative to each 
other. One can directly compare this information to the ac- 
tual positrons detected during SN 1987 A in order to assess the 
compatibility of the detectors, without ever referring to the ac- 
tual received spectrum of neutrinos. As explained below, we 
are assuming a common incoming neutrino spectrum. 

We summarize the visible energies (Ei) and experimen- 
tal uncertainties (A£y of positrons detected at the time of 
SN 1987A in Kam-II (11 data points marked with filled- 
diamonds) and 1MB (8 data point marked with open-circles) 
in the bottom panel of Fig. [TJ We exclude the Baksan Scin- 
tillator Telescope data 126l 12711 . which has much larger uncer- 
tainties associated with backgrounds, from our main results 
but comment on possible impact later. Data sets are synchro- 
nized based on the arrival time of the first events (note that 
we limit our study to time-integrated neutrino flux). While 
Kam-II originally reported 12 signal events, one of them, be- 
ing attributed to background, is excluded from our analysis, 
as in previous studies. The visible energy assigned to each de- 
tected event is based on the number of triggered photomulti- 
pliers from the Cerenkov light of the positrons. The statistical 
uncertainty on this visible energy is proportional to the square 
root of the total number of hits and is shown by a horizontal 
bar. 

Interestingly, the distributions of the detected positrons be- 
tween Kam-II and 1MB, as shown in the bottom panel, quan- 
titatively agrees quite well with our expectations based on the 
middle panel of Fig. Q] While Kam-II dominated the detec- 
tions below 16 MeV, the events between 16-30 MeV in each 
detector are comparable as one would deduce from the middle 
panel. At even higher energies (where the direct comparison is 
possible), 1MB detected 5 events while Kam-II detected only 
1 . If we accept that the 1MB observation of 5 events is a statis- 
tically probable outcome, since the effective number of targets 
in 1MB is almost twice to that of Kam-II at these energies, one 
would expect to see ~ 2.5 events in Kam-II. With this expec- 
tation, the probability of seeing < 2 events is over 50%, so that 
the observation of only a single event in the Kam-II sample at 
high energies is probably simply a downward fluctuation, and 
not in disagreement with the 1MB sample. Instead, if we as- 
sume that the Kam-II observation of 1 event is correct, then 
one would expect to see ~ 2 events in 1MB. With such a low 
expectation, seeing 5 events in 1MB is quite unlikely (with a 
15% probability of getting 4 or more events, and a 4% prob- 
ability of getting 5 or more events), suggesting instead that 
the true spectrum was indeed most faithfully sampled by the 
larger 1MB detector. 



Considering the low statistics and large experimental un- 
certainties involved, it seems that there is no obvious conflict 
between the distribution of detected positrons in the two ex- 
periments, suggesting the possibility of reconciling the Kam- 
II and 1MB data sets with a suitable incoming neutrino spec- 
trum. It is evident that while each experiment may suggest 
a radically different incoming neutrino spectrum individually, 
since the signal will be sampled distinctly by the detectors, 
only their proper combination could provide a sensitive probe 
over all energy ranges. In order to minimize the impact of 
statistics due to the small sample, we will present our subse- 
quent results with an emphasis on the combined data. 

IV. EXPECTED NEUTRINO SPECTRUM 

Neutrinos are produced thermally from the plasma and re- 
main trapped inside the collapsed core of a massive star un- 
til their optical depth becomes unity around the correspond- 
ing neutrino-sphere of each flavor. Then they stream through 
the envelope, preserving their initial energy distribution. The 
average energy at free streaming is dictated by the temper- 
ature at the neutrino-sphere, which is in turn determined by 
the strength of the coupling of neutrinos to the matter. One 
expects a hierarchy of energies, (E Vc ) < (Ep B ) < (E u „), 
since the electron neutrino flavor enjoys both charged and neu- 
tral current interactions (and couples to matter most strongly), 
while the non-electron flavors feel a weaker attachment to 
the plasma. The time-integrated luminosity and the hierar- 
chy of energies among neutrino flavors from a core-collapse 
supernova are still not well known and diverse values are re- 
ported by modelers |HL H 0, S 0, C3 Ell E3, E3l - While the 
supernova neutrino spectra are expected to be quasi-thermal, 
modifications due to non-standard effects, like neutrino mix- 
ing among various flavors J2IJ2I [30j pij |32j |33j [34l |35i 
neutrino dec ay |[38l[39ll40ll4lll42ll . neutrino-neutrino interac- 
tions additional channels of energy exchange 
between flavors ]5l|4J7j|4ja|49[], and/or any other novel mecha- 
nism due to unknown physics, may produce a time-integrated 
spectrum received on Earth that deviates significantly from a 
quasi-thermal shape. 

Many previous studies dealing with the SN 1987 A neutri- 
nos adopt a template neutrino spectrum and try to extract pa- 
rameters describing this spectrum from the Kam-II and 1MB 
data. Both Maxwell-Boltzmann 03|, <p(E u ) cx Ele~ E "/ T , 
and Fermi-Dirac spectra |51], (f(E v ) cx E^/[e E "^ T + 1], 
require only a single temperature parameter, T. An addi- 
tional degeneracy parameter, 77, is used to describe spec- 
tral pinching in a Fermi-Dirac distribution (52] . (p(E v ) cx 
E%/[e E "' T ~ v + 1]. More recently, a quasi-thermal distri- 
bution (based on a gamma distribution), is suggested 0], 
(p(E„) cx E2 e _ ( 7 + 1 )- E "/- E °, which can handle anti-pinched 
(broader) spectra, as well as pinched ones, through the pa- 
rameter 7. It has been shown that SN 1987 A data may favor 
a highly anti-pinched spectrum 115 311 . However, all such pre- 
defined functions describe only a very limited class of possi- 
ble shapes, which may not necessarily fit the actual received 
spectrum of neutrinos. 

More complex composite spectral shapes are also consid- 
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ered in the literature. Today, it is well established that neutri- 
nos may change flavor while they travel through matter, so- 
called MSW effect S HI H, H IM B- The den- 
sity profile of the supernova envelope, which neutrinos travel 
through, spans many orders of magnitude and thus enables 
many resonances and opportunities to swap flavors. The 1MB 
and Kam-II detectors are expected to observe almost identical 
spectra since the Earth effects l6ll l62l l63tl are unlikely to in- 
troduce significant differences (using recent determinations of 
the mixing parameters). A superposition of low-energy/high- 
luminosity and high-energy/low-luminosity thermal spectra 
due to neutrino mixing fits the data better than a single quasi- 
thermal spectrum Bill . A bimodal neutrino distribution the- 
oretically arising from the accretion and the cooling phases 
of a supernova B25I1 . is another example of a composite spec- 
trum. However, due to limited statistics, inferring multiple 
theoretical parameters is a challenging task, due to severe de- 
generacies among them. While both approaches effectively 
increase the number of parameters describing the spectrum, 
naturally improving the fit to the data, one can only probe an 
effective v e spectrum after any physical mechanisms modi- 
fying the spectrum, including oscillations, are taken into to 
account. 



V. NEUTRINOS FROM SN 1987A 

We resort to the data directly and attempt to find the sim- 
plest description of the effective v e spectrum that could be 



useful for studies which may require such input, regardless of 
the underlying physics and assumptions. As noted, we assume 
that the two detectors were exposed to the same incoming neu- 
trino spectrum. The theoretical relation between the received 
neutrino spectrum 4>(E U ) and the corresponding positron de- 
tection spectrum $(£' + ) is given in Eq. (0. The actual de- 
tected positron spectrum as recorded by an experiment can 
be expressed as a sum of "bumps" placed at the observations, 
Q(E + ) = J2i ~ Ei), where the index is over the set 

of events under consideration. This would be a faithful repre- 
sentation of the theoretically-expected positron spectrum (as- 
suming that the data was a probable outcome and not a statis- 
tical anomaly), only if the detected number of positrons were 
very large (i.e., in the case of a future Galactic supernova, op- 
erational detectors are expected to observe many thousands 
of events) and the uncertainties of the visible energies were 
small. Then, one can infer the received neutrino spectrum by 
inverting Eq. as 



4>(E V ) = 4ttD 2 J2 



5(E+ - Ej) 
N{Ei) a{E l + A) 



(3) 



which can be verified by substituting back into Eq. ©. When 
the statistics are limited, as here, we must regulate the dis- 
crete bumps in the spectrum by an appropriate smoothing; the 
most physically motivated method is to use the approximately 
Gaussian uncertainties on the measured energies. Hence we 
replace each S function by a Gaussian with a width of £j, 



5(E+ - Ek) - exp[-{E+ - Etf l{2e\ )]/(V2*e t ) . (4) 
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FIG. 2: The inferred neutrino emission spectra from either the Kam- 
II or 1MB data sets alone or their combination, as discussed in the 
text (taking into account the corresponding effective number of tar- 
gets and cross section). The shaded shape is a Fermi-Dirac spec- 
trum with canonical neutrino emission parameters (average energy 
Eo = 15 MeV and time-integrated luminosity L — 5 x 10 52 erg). 



Since very sharply peaked Gaussians would reveal the spuri- 
ous fine structure of the data, we choose a generous width of 
Ej = l.bAEi (where AEi is the uncertainty on the assigned 
energy of each detected positron), which enables us both to 
take into account uncertainties associated with the detection 
process adequately and to present a sufficiently smooth spec- 
trum. Any larger width would introduce excessive smoothing, 
which would obscure details of the distribution and spuriously 
enhance the tails. Note that in order to directly reconstruct the 
neutrino spectrum in this way, it is crucial that the neutrino 
and positron energies are related in a one-to-one way as in in- 
verse beta decay (unlike for neutrino-electron scattering, for 
example). 

Figure|2]displays the inferred spectra from either the Kam- 
II or 1MB data sets, as a function of neutrino energy. The 
differences between the inferred spectra from Kam-II or 1MB 
stems from a few facts. Since the 1MB detector has no sen- 
sitivity at lower energies, where we can rely on Kam-II only, 
the two least energetic events of 1MB, for which the efficiency 
is very low, then become disproportionately significant. Also 
the fact that only a single positron was detected by Kam-II at 
higher energies requires a significant suppression of the neu- 
trino spectrum. However, the effective number of targets for 
Kam-II is low compared to 1MB (as shown in the top panel 
of the Fig. [TJ and the expected number of events could easily 
fluctuate down. 

We have earlier established that Kam-II and 1MB probed 
different energy domains and also showed that the detected 
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positrons showed no obvious conflict with this expectation. 
So rather than focusing on the differences between the spectra 
suggested by the Kam-II or 1MB data sets, we will focus on 
what can be learned from the combined data set, since both 
experiments were measuring the identical incoming neutrino 
spectrum using an identical (water-Cerenkov) technique. 

By combining the two data sets and the corresponding (en- 
ergy dependent) effective number of targets (as displayed in 
the top panel of Fig. [TJ, we avoid over-emphasizing the dif- 
ferences due to low statistics. The combined result is also 
less susceptible to Poisson fluctuations associated with such 
limited statistics, yet covers the whole energy domain, and is 
more conservative. The neutrino spectra based on this com- 
bined data set is shown in Fig.[2](solid line). For comparison, 
we also show a Fermi-Dirac spectrum with canonical neutrino 
emission parameters (an average energy Eq = 15 MeV and a 
integrated luminosity L v = 5 x 10 52 erg). 

Since we construct the spectra directly from the data, it 
is not meaningful to make goodness-of-fit tests between the 
data and these constructed spectra (weighted with the energy- 
dependent proportionality factors). If we use delta functions 
to represent the data, then the cumulative distributions of 
the measured and 'predicted' spectra would be identical, so 
that a Kolmogorov-Smirnov test would indicate perfect agree- 
ment. When we use Gaussians to represent the data, the 
Kolmogorov-Smirnov test, while not strictly meaningful, pro- 
vides confirmation that our chosen width is not too large. 

While a pinched template spectrum puts more weight to the 
peak, and an anti-pinched spectrum puts more weight to the 
tail, the shape from the combined data can be best explained 
by a spectral shape that is enhanced both at the peak (to ac- 
commodate events from Kam-II where 1MB was not sensitive) 
and high energy tail of the spectrum (to better accommodate 
energetic events from 1MB), and depressed in between, com- 
pared to a thermal Fermi-Dirac spectrum. This basic shape 
of the underlying spectrum (see our Fig. 0, in agreement 
with the two-component composite spectrum of Lunardini and 
Smirnov could reconcile the Kam-II and 1MB data 

with each other. The corresponding luminosity for the com- 
bined spectrum, L v ~ 6 x 10 52 erg, is quite similar to that 
of the model, with an average energy, Eq ~ 12 MeV, that is 
somewhat lower. 

Since we consider a data set of 19 detected positrons, the 
overall uncertainty due to the Poisson nature of the detection 
will not be too large, l/\/l9 ~ 25%. In the peak region, 
our reconstruction is based on 12 events with a nominal un- 
certainty of 1/VT2 ~ 30%, so that the excess relative to the 
model is significant. In the tail region, there are 7 events with 
a nominal uncertainty of l/\/7 ~ 40%, which is also illus- 
trated by the differences between the spectra in this region of 
joint sensitivity. We emphasize that the tail region is not so 
uncertain, as there were 6 events above 35 MeV; this strongly 
precludes any hypothesized suppression of the tail. In Fig. [2l 
the combined spectrum above 40 MeV depends on the width 
chosen for the Gaussians, as there were no events at these en- 
ergies (however, there were 6 events between 35 and 40 MeV, 
each with relatively large energy uncertainties, and their sta- 
tistical weights must go somewhere). This is the least certain 
part of the spectrum, due to the low flux there. 



VI. DSNB DIRECTLY FROM SN 1987A 

Neutrinos from past core-collapse su pern ovae have been 
studied extensively IISHH MM, M, M, EzH IzE 
III H H [z! uSuEuEW and experimental limits 
are already suggesting an impending detection ll8ll l82l l83ll. 
The Diffuse Supernova Neutrino Background (DSNB) flux 
depends not only on the typical supernova neutrino emission 
spectrum, but also the expansion rate of the universe (redshift- 
time relation h(z) = yfJjvf^l + z) 3 + Ha, where =0.3, 
f^A = 0.7) and the core-collapse supernova rate (SNR), which 
presumably tracks the history of star formation rate (SFR) (see 
e.g. B84I1 and references therein). The neutrino emission per 
supernova is the most uncertain quantity 177118411 and is hence 
our focus here. The current precision of the data shows that 
the evolution of the SFR can be parametrized with a piecewise 
linear fit: 

Rsf(z) = R° SF {l + z) a for z<z p 

= TTT — J or z p< z < z ™<", (5) 

[L + Zp) 

where a ~ 3.44, (3 ~ 0, z p ~ 1 and z max ~ 5. 
The overall normalization in the local universe is R% F = 
O.OO95M /(yr Mpc 3 ). The core-collapse supernova rate as 
a function of redshift is R$n(z) = CHsf(z), while the frac- 
tion of stellar mass ending as supernovae, £ = 0.0132/Af Q 
(for the Baldry-GlazebrooklMF JUt]), depends on the stellar 
mass function. We note that this dependence on the assumed 
stellar mass function is small, as explained in Ref. [84]. While 
the directly measured core collapse supernova rate is slightly 
smaller [8(§[87l], recent studies suggest this may be mislead- 
ing. A large fraction of supernova exploding in very dusty 
starburst environments may go undetected, causing a 30-60% 
underestimate of the true core-collapse supernova rate B88I1 . 
Thus, while the supernova rate data are generally confirming, 
we use the more reliable star formation rate data for now. Ul- 
timately, it will be possible to use the supernova rate data to 
more directly predict the DSNB flux, as first pointed out by 
Ref. 117611 . The uncertainty on the flux due to astronomical 
inputs is small and will be further reduced with anticipated 
improvements in the data. In our analysis, we assumed that 
the normalization of the SFR and its evolution can be deter- 
mined independently by astronomers, and used as a fixed in- 
put to extract more accurate information on typical supernova 
properties. 

The neutrino emission per supernova times the supernova 
rate, when integrated over redshift and convolved with the 
cross section, yields the detected spectrum of DSNB neutri- 
nos, 

= ^{E v )N t d>(E v [l + z])^f^dz , 

Ho Jo K z ) 

(6) 

where c is the light speed and Ho is the Hubble constant. 
For a modern detector like Super-Kamiokande (SK), the ef- 
ficiency can safely be taken to be unity. We will present our 
results for a Super-Kamiokande sized detector of fiducial vol- 
ume 22.5 kton, corresponding to N t — 1.5 x 10 33 . 
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FIG. 3: The DSNB detection spectra based on the neutrino spec- 
tra inferred from either the Kam-II or 1MB data sets alone or their 
combination as in Fig. [2] compared to a model (shaded shape) with 
canonical neutrino emission parameters (the assumed core-collapse 
SN history is described in the text). 



In Fig. [3] the presented DSNB spectra are calculated by 
substituting the smooth neutrino spectra inferred from either 
the Kam-II or 1MB data sets alone or their combination as in 
Fig. |2 into Eq. ([§}. A canonical neutrino emission spectrum 
described by Fermi-Dirac distribution (E = 15 MeV and 
L = 5 x 10 52 ergs lT77to is also presented. 

Note that we have replaced each detected positron by a 
Gaussian to obtain a sufficiently smooth neutrino spectrum. 
This could inadvertently introduce some arbitrariness and dis- 
tort the results. Since the calculation of the DSNB already 
involves an integration over redshift, we can instead predict 
the DSNB directly from the data, which we call the (5-function 
prescription. When redshifting is taken into account, Eq. Q 
can be cast as: 



tf>(E v [l + z\) 



i 



6([l + z](E + +A)-(Ei + A)) 



5(z 



N(Ei) a(Ei + A) 
(Ei - E+)/E v ) 



E V N(E Z ) a(Ei + A) 



where we have used the delta function identity; 5 (ax — b) = 
S(x — b/a)/a. Substitution of this relation into Eq. (O elimi- 
nates the integration over redshift, so that the detected positron 
spectrum is 



V>(£+) = 



<j(E„)N t AirD 2 R SN ((Ei - E+)/E v ) 



H £f E lJ N(E l ) a(Ei + A) h((E t - E + )/E v ) 

(8) 

Figure |4] displays the resultant DSNB detection spectra, as 
deduced from the Kam-II or 1MB data sets alone or their com- 
bination, and compares to the model. Apart from the very high 



energy tail, both the smooth method (as presented in Fig. 0) 
and the 5-function method agree quite well over all energy 
ranges. In TableU we summarize the DSNB event rates in var- 
ious ranges of visible energy for both smooth and 5-function 
prescriptions (all values are quoted per 22.5 kton per year). 
The values in the last column, reporting the event rates in the 
energy range of 18-26 MeV can be compared to the event 
rate limit of ~2 at 95% C.L. (as inferred in Ref. 117711 from 
the v e flux limit of ~ 1.2 cm~ 2 s _1 at 95% C.L. reported by 
sk m through non-detection of excess counts above back- 
ground fluctuations), suggesting that the DSNB is already tan- 
talizingly close to detection. In the range 10-18 MeV, the 
event rates are nearly as large as the usual predictions based 
on supernova models, suggesting an imminent discovery of 
DSNB is well within the reach of current detectors, and espe- 
cially promising if thresholds are reduced by the addition of 
gadolinium II89I1 . 

There are various advantages to predicting the DSNB di- 
rectly from the data rather than through the intermediate 
stages of fitting functions. The DSNB mainly consists of the 
higher energy part of the neutrino spectrum and is much less 
susceptible to uncertainties/backgrounds around the threshold 
of detectors. While theoretically motivated formulas attempt 
to explain the bulk of the detected events, they tend to put 
more weight in the low energy part at the expense of intro- 
ducing distortions at the high energy parts of distributions 
relevant for the DSNB. For a high threshold of 18 MeV, al- 
most all of the statistical weight comes from the 1MB data. 
When the energy threshold is decreased, the Kam-II data will 
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FIG. 4: The DSNB detection spectra based on the neutrino spec- 
tra inferred from either the Kam-II or 1MB data sets alone or their 
combination, through the 5-function prescription which enables us to 
deduce the DSNB spectra without applying any smoothing (as dis- 
cussed in the text), compared to a model (shaded shape) with canon- 
ical neutrino emission parameters. 
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TABLE I: The DSNB event rates in various ranges of visible energy 
from the spectra displayed in Fig.|5]and Fig. [4] All quoted values are 
per 22.5 kton per year. 



Range (MeV) 




4-10 




10-18 




18-26 




S 


smooth 


S 


smooth 


S 


smooth 


Kam-II 


1.8 


1.7 


0.7 


1.0 


0.2 


0.3 


1MB 


0.6 


0.8 


1.4 


1.3 


0.6 


0.9 


Combined 


1.8 


1.8 


1.0 


1.2 


0.4 


0.5 


Model 




1.5 




1.8 




0.8 



be given more weight. Note also that we do not attempt to 
deconvolve any poorly-understood physical effects from the 
observed spectrum to deduce the original spectra at forma- 
tion. We are directly relating the observed spectrum of SN 
1987 A to the prediction for the observable DSNB spectrum 
from many supernovae. 

The idea of normalizing the DSNB flux prediction to the 
SN 1987 A data is appealingly empirical, and was first intro- 
duced by Fukugita and Kawasaki 02111 . and later developed 
in greater detail by Lunardini [22]. The disadvantage of re- 
lying on the SN 1987 A data is the uncertainties and appar- 
ent disagreements of the sparse data. While Fukugita and 
Kawasaki 112 ill used only the 1MB data, obtaining a neutrino 
emission per supernova in agreement with theoretical models, 
Lunardini showed that including the Kam-II data in the spec- 
tral fit leads to a more uncertain DSNB flux prediction, includ- 
ing the possibility of its being significantly lower than found 
by other authors. Our predictions for the combined Kam-II 
and 1MB data are consistent with those of Lunardini when 
differences in the star formation rate are taken into account. 
However, as we show, it is the 1MB data that are presently 
more relevant for the detectable DSNB flux, and accordingly, 
pessimistic DSNB detection rate predictions are disfavored. 
The fact that 6 of the SN 1987 A events were detected above 
35 MeV strongly precludes a supernova spectrum that is either 
too soft or too low. 

As with other forms of statistical inference, we are assum- 
ing that the observed data were representative of the truth (i.e., 
which is the maximally likely outcome), while considering 
how uncertain our subsequent results are. We are not using 
theory or other considerations to judge that the data were sub- 
ject to any particular statistical fluctuation. In particular, we 
found no evidence from the data alone that the Kam-II and 
1MB data sets were incompatible. A second point to con- 
sider is whether SN 1987 A was a representative supernova. 
Despite the large range of progenitors, and the wide vari- 
ety of optical supernova properties, the time-integrated neu- 
trino signals from core-collapse supernovae are expected to be 
mostly uniform, as they are well-connected to the properties 
of the newly-produced neutron stars. We note that both the 
SN 1987A data and the present SK limit on the DSNB flux 
both depend on statistical uncertainties at the ~ 30% level. 
It may well be that SN 1987 A is not a representative super- 
nova, and indeed, this is part of what we want to test with the 
DSNB 07711 . As shown in Figs. 3 and 4, and summarized in Ta- 
ble I, more typical results from theoretical expectations give a 



larger DSNB flux. Part of our point is that the SN 1987 A data, 
when considered without theoretical priors, do not support a 
very low DSNB flux. 

We now comment on the effects of considering also the data 
from the Baksan detector 1261 12711 . The relevant properties of 
this detector were very similar to that of Kam-II, except for 
being ~ 10 times smaller. Thus using the Kam-II yield, ~ 1 
event would be expected in Baksan, while 5 were observed, 
and with a somewhat higher average energy than in Kam-II. If 
the Baksan data were a faithful representation of the true spec- 
trum, then the corresponding DSNB flux would be enhanced 
by a factor > 5, which is likely already excluded. The most 
likely resolution is that Baksan saw an upward fluctuation of 
their large background rate, perhaps along with some events 
from SN 1987 A. 



VII. DISCUSSION AND CONCLUSIONS 

Neutrinos play a crucial role in the life and death of massive 
stars and so far, they are the only messengers that enable us 
to probe the inner workings of a core collapse. Thus, two 
decades after the detection of neutrinos from SN 1987 A, they 
continue to attract much attention. 

We assess the compatibility of the Kam-II and 1MB data 
sets in a model-independent way and turn to the data directly, 
skipping any intermediate stages and assumptions about the 
incoming neutrino spectrum. Our main conclusion is that if 
one drops theoretical prejudices on the spectral shape, the 
Kam-II and 1MB detections can be reconciled. The two data 
sets are primarily sensitive to different energy regimes of the 
neutrino spectrum (as seen in Fig. Q]), and only their proper 
combination probes the whole energy range. The actual ten- 
sion is between the adopted theory and the observations, since 
assumed theoretical shapes are always limited to a certain sub- 
set of all mathematical possibilities which do not necessarily 
describe the true underlying physics. When using the data 
directly, we consider the combined spectrum, based on both 
data sets, to be the most reliable. 

Using parameter-free inferential statistical methods, we 
have shown that the combined Kam-II and 1MB data can be 
best explained by a spectral shape that is enhanced both at 
the peak and the tail of the spectrum and depressed in be- 
tween, compared to a well-known Fermi-Dirac spectrum (e.g., 
Fig. [2]). Such a distribution may be arising from many differ- 
ent physical processes, including the details of core collapse 
or neutrino mixing, which are beyond the scope of this study. 
Once the effective received neutrino spectrum is determined, 
one can then work backwards and uncover the scenarios that 
will yield this measured spectrum. Considering that super- 
nova models still fail to robustly explode despite increased 
sophistication in modeling, it is an alluring possibility that a 
key element may still be missing 09011 . The model theoretical 
spectra with canonical emission parameters adopted in many 
studies could provide and capture essentials of the supernova 
neutrino spectrum, and be adequate for most purposes. While 
theoretical models still provide essential guidance, the neces- 
sity of fresh data on supernova neutrinos is obvious. The rarity 
of galactic supernovae, with the most optimistic rate estimates 
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of at most a few per century, makes this a challenge. 

Apart from the proposal to detect individual neutrinos from 
galaxies within the 10 Mpc neighborhood of the Milky Way 
with future Mton-scale detectors 19111 . the DSNB presents 
the only sensible alternative to gain information on neutrino 
emission from supernovae. We use a nonparametric approach 
in order to determine one observable, the Diffuse Supernova 
Neutrino Background, directly from another observable, the 
SN 1987A data on supernova neutrino emission, rather then 
proceeding through intermediate stages of fitting functions. 
We show that this prediction cannot be too small (especially 
in the 10-18 MeV range), since the majority of the detected 
events from SN 1987 A were above 18 MeV (with 6 above 
35 MeV). We emphasize that our DSNB predictions are not 
very dependent on the details of our analysis procedure, as 
can be seen by comparing Figs. [3] and |4] and especially by 
comparing the "6" and "smooth" cases in Table I. 

These results are nearly as large as the usual predictions 
based on supernova models, suggesting an imminent discov- 
ery of DSNB is well within the reach of current detectors, 



and especially promising if thresholds are reduced by the ad- 
dition of gadolinium M89I1 . A gadolinium-enhanced Super- 
Kamiokande should also be able to provide a measurement 
of the spectral shape lf77l l92ll . providing further clues. The 
DSNB, which may even one day serve as a tool for extracting 
cosmological evolution parameters [93], is a leading candi- 
date along with the other contenders, like cosmogenic neutri- 
nos 09411 . to open new doors to the cosmos and provide the first 
neutrino detection originating from cosmological distances. 
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